Equivalence relations for the classical capacity of single-mode Gaussian quantum channels 
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We prove the equivalence of an arbitrary single-mode Gaussian quantum channel and a newly Ae,fme,A fiducial 
channel preceded by a phase shift and followed by a Gaussian unitary operation. This equivalence implies that 
the classical capacity and Gaussian classical capacity of any single-mode Gaussian channel can be calculated 
based on this fiducial channel, which is furthermore simply realizable with a beam-splitter, two identical single- 
mode squeezers, and a two-mode squeezer. In a large domain of parameters, we also provide an analytical 
expression for the Gaussian classical capacity, exploiting its additivity, and prove that the classical capacity 
does not exceed it by more than 1/ In 2 bits. 



Introduction. — Quantum channels play a key role in quan- 
tum information theory. In particular, bosonic Gaussian chan- 
nels model most physical communication links, such as opti- 
cal fibers or free space information transmission |TJ[2]]. One 
of the central characteristics of quantum channels is their clas- 
sical capacity, i.e. the maximal amount of classical bits that 
can be transmitted per use of the channel. A lot of attention 
has already been devoted to the study of the classical capac- 
ity of Gaussian channels [3— 23 ]. Since Gaussian encodings 
are more relevant for experimental implementations, easier to 
work with analytically, and conjectured to be optimal lfl6l . 
the so-called Gaussian classical capacity was introduced and 
evaluated for specific Gaussian channels lfl6H23l . 

In this Letter, we greatly simplify the calculation of these 
capacities [42 1 for an arbitrary single-mode Gaussian chan- 
nel. Namely, we show that any single-mode Gaussian chan- 
nel is indistinguishable from a newly definedfiducial channel, 
preceded by a phase-shift and followed by a general Gaussian 
unitary. Since neither the phase-shift at the channel's input nor 
the Gaussian unitary at the channel's output change the output 
entropy or affect the input energy constraint, we conclude that 
the capacities of this channel are equal to those of the fidu- 
cial channel, which furthermore, is completely specified by 
three parameters. This conclusion also holds for any cascade 
of Gaussian channels since the latter is equivalent to another 
Gaussian channel. Our results allow us to go beyond previous 
works on the Gaussian capacity [20-22) and provide its an- 
alytical expression for any Gaussian channel in some energy 
range, where it is additive. Then, we prove that the classical 
capacity does not exceed the Gaussian capacity by more than 
l/ln2 bits (generalizing )24|), which furthermore becomes 
the actual capacity if the minimum output entropy conjecture 
for phase-insensitive Gaussian channels |fT6ll25l is true. 

Gaussian channel. — Let p G (a,V) be an n-mode Gaus- 
sian state, where the coherent vector a. 6 M. 2n and the 
2n x 2n covariance matrix (CM) V are the first- and second- 
order moments of the 2n dimensionless quadratures, respec- 
tively, with h = 1. Then, a Gaussian channel <!> is a 
completely-positive trace-preserving map which is closed on 
the set of Gaussian states [16]. It transforms Gaussians states 
as $[/3 G (a in , V in )] = p G (a out ,V out ), with the input mo- 



ments {c*in,Vin} being connected to the output moments 
{a out , Vout} according to 

a out = Xa in + S, V out = XV in X T + Y, (1) 

where 8 is the displacement introduced by the channel, X is 
a 2n x 2n real matrix, and Y is a 2n x 2n real, symmetric, 
and non-negative matrix. For simplicity, we choose 6 = in 
what follows (the capacity is not affected by first-order mo- 
ments), and focus on the action of the map $ on second-order 
moments using the simplified notation ^(Mn) = Vout- Then, 
the map $ is fully characterized by matrices X and Y, which 
must satisfy the relation Y + § (ft - Xf!X T ) > |fl6l , 
where 

°-©(-°ii) < 2 > 

is the symplectic form |26l . In the following, we focus on 
single-mode channels (n = 1) and use the parameters 

t = det X, y = VdetY, (3) 

where r may be a channel transmissivity (if r < 1) or ampli- 
fication gain (if r > 1), while y reflects the added noise. The 
map $ describes a quantum channel if y > ||r — 1| 1271 . and 
it is an entanglement breaking channel if y > |(|t| + 1) l28l . 
The single-mode Gaussian channels can therefore be conve- 
niently represented in a (r,y) plane, see Fig.[T] 

Canonical decomposition. — Any single-mode Gaussian 
channel $ can be decomposed as $ = U 2 ° <1> C o U\, where 
U\ and U2 are Gaussian unitaries, and $ c is a canonical chan- 
nel characterized by the matrices (Xq, Yq) [29-31]. The ac- 
tion of a Gaussian unitary U on a Gaussian state can be com- 
pletely specified by a symplectic transformation M acting on 
the second-order moments of the state (we ignore first-order 
moments), so that the canonical decomposition may be writ- 
ten as (U 2 o $ c o l7i)(V5n) = M 2 $ c (MiM„M 1 T ) Mj. It 
is known that there are seven classes of canonical channels 
$ c (see Table |l| 129H3T1 . We first discuss these canonical 
channels $ c and give explicit relations for the canonical de- 
composition of an arbitrary Gaussian channel <i> (see also the 
supplemental material [43 1). 
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FIG. 1. Admissible regions in the parameter space (r, y) for Gaus- 
sian quantum channels. Each thermal channel &™ y ) is associated 
with a point. The vertical line at r = corresponds to the zero- 
transmission channel as well as the classical signal channel 4> cs . 
The vertical line at r = 1 corresponds to the classical additive-noise 
channel. The limiting case (t = 1, y = 0) is the perfect transmission 
channel as well as the single quadrature classical noise channel <J> SCJ . 
The Gaussian capacity of 5>f T y s \ is additive if TV > TVthr. This is 
equivalent to y < y thl = |rj(e- 2s (l + 27V) - 1)/(1 - e~ is ). An 
example of j/ t hr is given by the dashed line, for particular values of 
N and s. 



The first five channels in Table|I]can be treated together, and 
we refer to them collectively as thermal channels, denoted by 



$7?..y They are defined as 



X 



TH 









sgn(r) 



'- TH 



y o 

y 



(4) 



where sgn(r) = — 1 if r < and sgn(r) = 1 if r > 0. 
As shown in Fig. [2] (a), the channel $ TH can be physically 
realized by a beamsplitter with transmissivity T followed by 
a phase-insensitive two-mode squeezer (TMS) with gain G 
0271 . For the zero-transmission (r = 0), lossy (0 < r < 1), 
amplification (r > 1), and classical additive-noise channel 
(t = 1), the output is given by the signal's output of the 
TMS, and these four canonical channels in Table [^correspond 
to phase-insensitive channels. For the fifth canonical chan- 
nel, i.e., the phase-conjugating channel (r < 0), the output 
is given by the idler's output of the TMS. These five chan- 
nels share the property that they map any thermal state to a 
thermal state, therefore, we call them thermal channels. Each 
particular channel $ TH with parameters (r, y) corresponds to 
a single point in Fig.[T] while the relations between (t, y) and 
(T, G) are given in Table [i] 

The sixth canonical channel is the classical signal (or 
quadrature erasing) channel, which we denote by $ cs . Its 
action is defined by 
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This channel can be physically implemented with a conti- 
nuous-variable controlled-NOT (C-NOT) gate l32ll33l . The 
corresponding scheme is depicted in Fig. [2] (b), where 
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FIG. 2. Physical realization of (a) the thermal channel <E> , (b) 
the classical signal channel $ cs , (c) the single-quadrature classical 
noise channel $ SCJ , and (d) the fiducial channel $ F . T stands for a 
beamsplitter with transmissivity T, G for a two-mode squeezer with 
gain G, S for a single-mode squeezer, j 0) for the vacuum state, and 
"H" denotes "tracing out" the mode. 



G = y + 1/2. Note that r = 0, implying that $ cs is al- 
ways entanglement breaking (see Fig. [TJ. The seventh and 
last canonical channel is the single-quadrature classical noise 
channel, which we name $ S Q. Its action is described by the 
matrices 



^SQ — 




\ 



(6) 



which corresponds to r = 1 and y = (as the perfect trans- 
mission channel). It is known that <1> S Q is not entanglement 
breaking 11281 . It is shown in supplemental material ll43l that 
channels $ cs and <£> S Q can be represented by schemes Fig. [2] 
(b) and (c), respectively. 

Fiducial channel. — Now, our central point is that the above 
canonical decomposition is not always useful for evaluating 
capacities of bosonic channels with input energy constraint 
(which is needed, otherwise the capacities are infinite). In- 
deed, the Gaussian unitary U\ that precedes the canonical 
channel <£> c changes, in general, the input energy. There- 
fore, we introduce a new decomposition in terms of a fiducial 
channel <£> F , where the preceding unitary is passive and does 
not affect the input energy restriction. Below we show that 
this decomposition has the major advantage that the energy- 
restricted capacity of any Gaussian channel reduces to that of 
the fiducial channel $ F . The latter generalizes <£> TH by adding 
a squeezing to the effective added noise, that is 







"-TH, 



-2s 



(7) 



Thus, it depends on three parameters (r, y, s), and will be de- 
noted by $ F T y g » . This channel can be physically realized 
by the setup depicted in Fig. [2] (d), where the "idler" corre- 
sponds again to the output of the phase-conjugating channel 
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TABLE I. Canonical channels $ c with corresponding symbols as defined in 129443 1 1 . new representation in terms of $ TH , $ cs , <J? SCJ and 
the corresponding matrices (Xc, Yc), where a z = diag(l, —1). The parameter T £ [0, 1] is the transmissivity of the beamsplitter and the 
parameter G > 1 is the gain of the two-mode squeezer shown in the physical schemes in Fig. [2] 



and the "signal" to that of the other channels. In the case of 
< t < 1, this channel corresponds to the mixing of the in- 
put state with an arbitrary squeezed thermal state on a beam- 
splitter with transmissivity r. The fiducial channel <D F can 
be used to decompose any Gaussian channel $ (at least in a 
proper limit): 

Theorem 1. For a single-mode Gaussian channel $ defined 
by matrices X and Y with t / and y > 0, there exists a 
fiducial channel $ F defined by matrices Xp(r), Yp(y, s) with 
t and y obtained from Eq. p|, a symplectic transformation 
M, and a rotation in phase space such that 

X = M X f (t)&, Y = MY F (y,s)M T , (8) 

where the explicit dependencies of M, ®, and s on the pa- 
rameters of the channel $ are presented in the supplemental 
material in Eqs. $26) - $56) . 

Proof. We only sketch the proof here (see ll43ll for the 
full proof). First, one finds matrices ®y and Sy such 
that S^QyYQySy 1 = diag(y,y), where and S de- 
note matrices corresponding to a rotation and a symplec- 
tic squeezing operation, respectively. Secondly, one obtains 
the singular value decomposition X = ®i X SxXp®2x, 
where Xp reads as in Eq. (|7j. Then one defines M = 
®i X S x ®l, where F is found such that M l YM = 
Yp = y diag(e 2s , e~ 2s ). The squeezing parameter s depends 
on all angles and squeezing operations Sx , Sy ■ Finally, one 
introduces 0p in X, i.e. X = ®i X Sx®f®fXf®2x = 
MXp®, where © depends on ®i X ,®f and the sign of 
r. Despite Theorem [T] requires that r =^ and y > 0, or, 
equivalently, that rank(X) = rank(l^) = 2, it can be ex- 
tended to lower-rank cases with minor modifications as shown 
in El. □ 

Capacities. — We now focus on the capacity C of the Gaus- 
sian channel $ with a constrained mean photon number N at 
the input. It is defined as the maximal amount of bits that can 
be transmitted per use of the channel $, i.e. |[T5ll34l 

C($,iV)= lim -C x ($^™,nA>), (9) 



where n is the number of channel uses, and C x ($, N) is the 
one-shot capacity of the channel, i.e. 

C x ($,A0= max x(*,aO, 

fj, : pGf jv 

r (10) 

x(*,l*) = S(p\p\)- J li(dx)S{*\p x ]). 

Here S(p) — — Tr(/51og 2 p) is the von Neumann entropy. The 
maximum is taken over all probability measures li(x) in the 
whole space H of pure states p x of the field mode and p — 
J p(dx)p x belongs to the set £ x of states which have a mean 
photon number not greater than N. Since it was shown that, in 
general, the one-shot capacity is not additive [35], one has to 
take the limit in |9]), unless additivity is explicitly proven for 
the given channel. The decomposition stated in Theorem [T] 
implies: 

Corollary 1. For a single-mode Gaussian channel $ with pa- 
rameters (t ^ 0,y > 0), there exists a fiducial channel $ F as 
defined in Theorem^such that 

C{$,N) = C($ F ,7V). (11) 

Proof. First, the symplectic transformation M that follows 
<£> F in Theorem [I] does not change the entropies in the def- 
inition of x- Second, there is no energy constraint on the 
output of the channel. Hence, M can be dropped. Then, 
the rotation preceding $ F in Theorem fT] may be re- 
garded as change of a reference phase that can be chosen 
arbitrarily; therefore, can be dropped as well. Thus, 
C x ($, N) = C X ($ F , N) holds. To evaluate the one-shot ca- 
pacity of n copies of $ the same reasoning can be applied, 
where now the preceding and following transformations are 
given by ®" =1 M and ®f =1 0, respectively. Hence, it follows 
that C x ($® n , nN) = C X (($ F )®", nN) which together with 
Eq. |9]) implies Eq. ( fTT| . Note that despite Eq. ( fTT) requires 
t / and y > 0, it can be easily extended to the general case 

EH. □ 

Corollary [T] has also other implications. If the correspond- 
ing fiducial channel $ F r y ^ is entanglement breaking, i.e. 
U > (l r I + l)/2, then its one-shot capacities of both <£> F T y ^ 
and $ are additive I1361I371 . namely C($, N) = C X ($ F , N). 
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Gaussian capacities. — For experimental implementations 
and analytical calculations, it is useful to focus on Gaussian 
symbols and encodings. We call the capacity restricted to 
Gaussian encodings the Gaussian capacity [20-22] C G and 
define it as 



C G ($, AO = lim ^C G ($ 



n— foo Tl 

max 



nN), 
X($,M G ), 



(12) 



where C G ($,N) is the one-shot Gaussian capacity. The 
maximum is taken now over all measures /i G (a, V) for which 
p G (a in , V in ) = J p G {da, dV)p G (a, V) is in the set of 
Gaussian states with a mean photon number not greater than 
N. Unlike previous works (e.g., ifTTlO . we require the individ- 
ual symbol states as well as the averaged state to be Gaussian. 
Then we prove that the one-shot Gaussian capacity of an ar- 
bitrary single-mode Gaussian channel $ is given by the well 
known expression |4|: 



C G (*,N)= max {f{v,v) I Tr[M„ 



V mod ] < 2A>+1}, 



x G (v,v) = g 



(x + 1) \og 2 (x + 1) - x log 2 x, 



where VJ n is the CM of a pure Gaussian input state p G (0, V n ) 
fulfilling det (2Vi n ) = 1, Vnod is the CM of a classical 
Gaussian distribution used to displace the input state and to 
generate the modulated input state p G (0, V n ) whose CM is 
V in = V hl + Vnod with Tr[V m ] < 2N + 1. Furthermore, 
v = \/det V ou t and v — \/ det V^ ut are the symplectic eigen- 
values of the output and modulated output states with CM 
V ont = <£(Vin) and V out = $(V in ), respectively. The de- 
tailed derivation is provided in |43 1. 

The one-shot Gaussian capacity is equal to the Gaussian 
capacity, i.e. C G (<1>, N) = C G ($, N), if it is additive. Inter- 
estingly, such an additivity can be proven if the input energy 
exceeds some threshold A^r (see (43j). Note that ifTTll also 
derives additivity but for a slightly different definition of C G 
and without respecting the energy threshold. In addition, an 
analog of Corollary [T] can easily be shown to hold for Gaus- 
sian capacities, namely C G (<1>, N) = C G ($ F ,N). There- 
fore, using the fiducial channel <£> F , we can analytically find 
the Gaussian capacity of any Gaussian channel in this high- 
energy regime: 

Corollary 2. For a single-mode Gaussian channel <f> with pa- 
rameters (t ^ 0,y > 0), there exists a fiducial channel $ F as 
defined in Theorem^such that 

C G ($,N) = C G (<Z>J TtViS) ,N) 



g ( \t\N + y cosh(2s) + 



1 



T 1 



g\y + 



T 1 



ii' V ■ X llu = | [ e 2 l s l + ^sinh(2| S |)-l 



(14) 



The proof is presented in [43 1. Note that the energy thresh- 
old TVthr depends on the parameter s characterizing the fidu- 
cial channel $ F v. For thermal channels <I> TH = $ F m , 
the threshold N tnr vanishes and the additivity holds on the 
entire energy range. Then, Eq. ( fT4| ) coincides with previously 
derived expressions for particular cases l4l [T6ll . In Fig. [T] we 
illustrate an example of the domain where A~ > Athr> hence 
where Eq. < fl4] > holds. Note, that Eq. ( |T4] > becomes the actual 
capacity C($, N) (for A" > Athr) of an arbitrary single-mode 
Gaussian channel $ provided that the vacuum state is proven 
to minimize the output entropy of a single use of an ideal am- 
plification channel l27l[38l . 

Upper bounds. — Recently, upper bounds have been derived 
on the capacity of phase-insensitive channels, i.e. <I> TH with 
r > l24ll39l . Using Corollary [2] we can generalize it to any 
Gaussian channel in this high-energy regime: 

Corollary 3. For a single-mode Gaussian channel $ with pa- 
rameters (r > 0, y > 0) and N > A^hr, 



C G ($, N) < C($, N) <C < C G ($, N) + 



ln2' 



/2rA^+(22/ + l-r)sinh 2 s 
y V 2y + 1 + t 



(15) 



(^3) where C G ($, N) is stated in Eq. fB] ). 



Proof. The fiducial channel corresponding to <f> can by defini- 
tion be decomposed into $f TiJ/iS ) = g=i s) ° $ F T i^r s} 
with T = 2r/(2y + r + 1) [see Fig.[5]and Table|I|. Then, we 
can upper bound the capacity of <£> F r ^ by the capacity of the 
first channel, i.e. 



*^ sV N)<C, 



where C = g(TN + (1 - T) sinh 2 s) ED. Then we define 
A(s) = C-C G = g[A{B+l)- l ]-g{A+B cosh 2 s)+g(B), 



where A = tN + [y - ^=1] sinh 2 s and B = y + T -^. It 
was shown in J24) that A(0) < l/ln2. Since A(s) < A(0), 
Vs, the Corollary is proven. □ 

Note that for r < we can state a similar upper bound on 
the capacity, C($, A^) < C, where C is given by Eq. ( p~5] > 
with the replacement y —> —y. However, in this case the last 
inequality in Eq. ( fT3] l does not hold. The upper bounds of 11391 
can also be extended in a similar fashion as shown in [43 1. 

Conclusions. — We have shown that an arbitrary single- 
mode Gaussian channel is either equivalent to a newly de- 
fined fiducial channel preceded by a phase-shift and followed 
by a Gaussian unitary, or can be obtained in a proper limit 
of this combination. This is exploited to reduce the energy- 
constrained classical capacity of any single-mode Gaussian 
channel to that of the fiducial channel. We give an analytical 
solution to the Gaussian capacity above some energy thresh- 
old, where additivity can be proven, and show that in this case 
the classical capacity does not exceed it by more than l/ln2. 
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Supplemental Material 

Physical representation of channels <J? CS and $ SQ 

In the following we explain the physical schemes of the channels $ cs and $ s< ^ as depicted in Fig. [5] (b) and (c). First, 
we discuss the main "building block" of these schemes, namely the CV-CNOT gate acting on a two-mode state with CM 
Vi n th = Vin © Kh> consisting of a general input mode with CM V sn and an ancillary mode being in a thermal state with CM 

v m =( v « v A, v t ^d °y (i6) 

\Vqp v p J \0 yj 

The ancillary mode with CM Vth in case of the channel $ cs is in fact the output of the two-mode squeezer with gain G [Fig. [2] 
(b)], after having traced over the second ancillary mode, therefore, y = G — 1/2. The action of the symplectic transformation 
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FIG. 3. Equivalence of (a) an arbitrary Gaussian channel $, (b) the canonical decomposition containing a canonical channel $ c and (c) the 
decomposition in terms of the fiducial channel $ F as stated in Theorem[T] 



S'cnot corresponding to the CV-CNOT on the joint state reads 
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Tracing over the first mode [as shown in Fig. [2] (b)] leads to the output for the classical signal channel $ cs . Tracing over the 
second mode [as shown in Fig.|2](c)] and taking the ancillary mode to be in a pure state (y ~ 1/2) leads to the output of the 
single quadrature additive noise channel $ S Q. 



Proof of Theorem [T] and method how to obtain new decomposition 

Theorem 1. For a single-mode Gaussian channel $ defined by matrices X and Y with t =/= and y > 0, there exists a fiducial 
channel <!> F defined by matrices Xp(r), Yp(y, s) with r and y obtained from Eq. a symplectic transformation M, and a 
rotation in phase space © such that 

X = MX f (t)&, Y = MY F (y,s)M T , (18) 
where the explicit dependencies of M, 0, and s on the parameters of the channel $ are presented in £<7^.( |26| >-( |36| >. 
Proof. The action of the channel $ on an input CM V reads as in Eq. ([TJ 

<I>(V) = XVX J + Y, (19) 
where X is a real 2x2 matrix and Y a real, symmetric and non-negative 2x2 matrix. In the main text we defined the parameters 

t = det X, y = Vdet Y, (20) 

which have to satisfy 

II 



> 



(21) 



in order for the map to be physical. In [29-31 1 it was stated that for any Gaussian channel $ there exists a canonical decompo- 
sition U2 ° ( E >C o Ui, where $ c is a map belonging to one of the seven canonical types that are stated in the main text in Table|l] 
The corresponding action on the CM reads 



*(V) = M 2 (X c M 1 VMjX c + Y c )Mj, 



(22) 



where Xq, Yq are the matrices defining the canonical channels (see Table [I] in the main text). In the following we obtain the 
new decomposition in terms of the fiducial channel as stated in the Theorem and furthermore, confirm Eq. ( p2] >. The proof is 
structured as follows. For given matrices X, Y we have to distinguish three cases which depend on the ranks of X and Y and 
correspond to canonical decompositions for which $ c is either $ TH , <1> S Q or $ . In the first case our new decomposition will 
contain finite squeezing operations, while for the other two cases the new decomposition is shown to be valid in a proper limit 
of infinite squeezing. 
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For the following calculations we define the symplectic matrices corresponding to a rotation and a squeezing operations, i.e. 

o(U) = r ""' )• ^> = i • °V (23) 




where we sometimes omit the explicit dependence on the rotation angle or squeezing parameter. For the given CM Y one can 
find a rotation 0y = O(0y), such that 0}Y0y = diag(j/i, y 2 ), where y ll y 2 > are the eigenvalues of Y. Since matrix X 
is always real it has a singular value decomposition (SVD) 

X = ® 1X A X J@ 2X , (24) 
where ®i X = 0(6\x) an d ®2x = O(0 2 x) are rotation matrices. Here 

a , • /■ ^ r / 1 ifr>0 

A x = diag(xi,x 2 ), J=< , (25) 

\ a z if r < 



where x\,x% > are the singular values and a z — diag(l, —1). Using equality detX = det (Ax J) and Eq. p0[ ) we get 
t = ±xiX2 and y = -Jyyyi. The condition on the determinants of X and V stated in Eq. pTj ) allows us to exclude the 
following combinations of ranks because they are non-physical: (rank(X), rank(l r )) ^ {(0, 0), (0, 1), (1, 0), (1, 1)}. The 
physically allowed combinations of ranks therefore read (rank(X), rank(l^)) € {(0, 2), (2, 0), (2, 2), (1, 2), (2, 1)}. Below we 
treat the first three "physical" couples together and the other two individually. 

1) First, we derive the relations for the case rank(X) = rank(Y") = 2. The latter implies that X\, x 2 , yi, J/a 0. Then we 
can construct the squeezing operation Sy = S(sy), with sy = \ In (2/1/2/2) such that diag(?/i, y^Sy 1 — diag(y, y). This 
implies that 

F = Y S y r TH Sy0 y = y®YS Y &y, Y TH = diag{y,y). (26) 

Here the symplectic transformation &ySy realizes the symplectic diagonalization of Y, where y is the symplectic eigenvalue. 
Furthermore, we can define a squeezing operation Sx = S(s x ), with sx = \ In {xx/x 2 ), such that Eq. ( p4| ) can be written as 



X — ® lx S x X Tii ® 2x , -^th — I v „ . 1 /;— r] • ( 27 ) 

^ sgn(r)V|r|y 

Notice that the matrix Xth has the property that 

X TH O(0) = O(sgn(r)0)X T H. (28) 

Now we multiply Y in Eq. |26) from both sides with the identity matrix in the form I = M® F S X ©J x , 

IYI = yM& F S x 1 ®l x &YS Y &l& 1 xS x 1 @lM T , (29) 

where 

M = ® 1X S X ®1. (30) 
Here 0p is a rotation that must be chosen in way such that matrix Yp is diagonal, i.e. 

Y F = y® F S x 1 ® T lx ® Y Sy®l® lx S x 1 ®l = j/diag(e 2s ,e- 2s ), (31) 
which leads to Y = MY F M T . The explicit expression for the squeezing parameter s reads 



£ = (1 + e- 4s n(l + e 4sx ) + (-1 + e- 4 ^)(-l + e 4sx ) cos(2(# y - 9 1X )), 
£ = (1 + e- iSY ){-\ + e isx ) + (-1 + e- 4 ^)(l + e 4sx ) cos(2(#y - 1JC )). 

The angle 6>f of rotation 0p = 0(9 F ) reads 



(32) 



arcsin f jp£L ) , (33) 
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where 

= e~ 2 ^(|- y/-16e^- ay ) + g2) 

2sin(2(0y-0 lx ))(-l + e-^) ' ^ ; 

Using definition Xp = Xth [see Eq. |7| in the main text] and Eq. ( p8| one can rewrite Eq. p7| ) as 

X = © 1x S x ©t 0fXf 2x = lx 5 x ©;X F O(sgn(r)0 F )© 2x = MX F 0, (35) 

where 

= O(sgn(r)0 F + 2x ). (36) 



In summary, we found that 



X = MX F e, Y = MY F M T , (37) 



and we have proven the theorem. 

Now we extend the theorem to other combinations of ranks of X and Y, Let us consider the case rank(X) = 2, rank(Y) = 
0. Since Y = it follows that y = 0, which together with Eq. pi) implies that r = 1. Therefore, the channel is unitarily 
equivalent to the perfect transmission channel. All relations derived above are found in the same way where one has to fix 
sy = 9y = 0, which leads to Sy = 0y = I. 

The third case corresponds to rank(X) = 0, rank(Y') = 2. It can also be treated using the above relations. In this case 
X = 0, r = which together with Eq. pi) implies that y > 1/2. This channel is unitarily equivalent to the zero-transmission 
channel and has trivially a capacity equal to zero. The decomposition containing the fiducial channel is found above where one 
has to fix sx — dix — @2x = 0. 

We remark that by inserting Eqs. |26| and ( [27] ) in Eq. (JT9J one obtains the canonical decomposition $ = [7 2 o <1> TH o U\, 
which in terms of the symplectic transformations reads as in Eq. p2|), with 



X c =X T h, Y c = Y tu , M 1 = Sy 1 @' Y T ® lx 'Sx®2x, M 2 = & Y S Y , (38) 

where 0^ = 0(sgn(T)6> y ) and ® lx ' = 0(sgn(r)6»i X ). This confirms that for (rank(X), rank(^)) € {(2, 2), (0, 2), (2, 0)} 
the physical action of $ indeed corresponds (up to unitaries) to the action of $ TH . In Fig. [3] we sketched the equivalences found 
above. 

2) Next, we consider the case rank(X) = 2,rank(l^) = 1. This implies y = and together with Eq. ( |2~Tj i that r = 1. 
The eigenvalues of Y now read y\ — 0, 2/2 > (the case y% > 0, y 2 = follows the same treatment). Similarly to the case 
rank(i') = 2 one can find a rotation 0y such that ® y Y®y = diag(0, y 2 )- Then, one can construct a squeezing operation 
®y with sy — — \ ln(2y 2 ) which yields 



Y = ® y S y Y sq S y ®y, F S Q = diag^O,0 ; 



(39) 



The matrix Y$q can be recovered with an additional squeezer St — S(st) in the limit of infinite squeezing, i.e. Y$q = 



liirisy^oo he 2st S t 2 from which follows 

Y= lim \e- 2sT ® Y S Y S T '' z SY®Y = I™ \e- 2sT ® Y S YT ® Y , (40) 

where Syt = S(sy — st)- Since rank(X) = 2 we can decompose X as in Eq. < |27| > but with the simplification r = 1, i.e. 

X = ® 1X S X X SQ ® 2X , X sq =L (41) 

We observe that we can replace X^q = Xp, where X F is as defined as above with r = 1. Thus, we get the same decomposition 
as stated in Eq. ( f27j ). Now one can recover both matrices X, Y as a limiting case of Eq. ( [37] ), namely, 

X = lim MX F 0, Y = lim MY F M T , (42) 

where in the definitions of iVf ( |30[ l, ([36) and Yp ([3TJ in one has to make replacements r —> 1, sy — > s Y — st and y —> \e~ 2sT . 
We remark that this replacement only effects matrix Yp and rotations 0p and 0. Thus, we found both matrices X and Tasa 
limiting as of the decomposition stated in the theorem. 
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Note that by inserting Eqs. ( |4T} and ( [39} into Eq. ( fT9} we find the canonical decomposition $ = [/jo <I> S Q o f/i, which in 
terms of the symplectic transformations is given by Eq. d22|, with 



X C =X SQ , Y C =Y SQ , M 1 = Sy 1 ®i® lx S x @ 2x , M 2 = @ Y S Y , (43) 

showing that the physical action of $ in this case corresponds (up to unitaries) to the action of $ S Q. 

3) The last case to consider is rank(X) = 1, rank(Y) = 2. Since in this case r = it follows from Eq. pT} that y > i. 
The SVD of X now reads X = ©ixdiag(xi, 0) (the case X\ — 0,^2 > follows the same treatment). One can define 
Sx = S(sx) with sx = ln(a;i) such that 

X = & 1X S X X CS , X cs = diag(l, 0). (44) 
Since Xcs can be expressed as Xcs = lini ST ^oo e~ ST St, where St — S(st), Eq. ( |44-[ > becomes 

X = lim e- ST & lx SxS T = Km e~ ST & lx SxT, (45) 

where Sxt = S{sx + «t)- Since y > \ we find as in the case y > treated above [see derivation of Eq. |26|], a rotation ® Y 
and squeezing SV such that 

Y = & y S y YcsSy®y, Y cs =di a g(y,y), y> 1 -. (46) 
We observe, that we can recover matrices X,Tasa limiting case of Eq. p7] i, i.e. 

X = lim MX] €). Y = lim MY F M T , (47) 

st— ^oo st— >oo 

where in the definitions of M ( [30] > and ( |36| ) one has to make replacements Q 2 x ~> 0, sj — > + and r — >• e _2sT . Note 
that this replacement effects Af but does not effect matrix Y" stated in Eq. ( |46] i. Therefore, we found also for the last case both 
matrices X, Y as limiting cases of the decomposition stated in the theorem. 

The canonical decomposition in this case is found as follows. By inserting Eqs. ( |46| i and (|44j> in Eq. ( [T9] >, we obtain 

*(V) = M(XVX T + Ycs)M T , X = S^SyQixSxXcs, M = ® y S y . (48) 

For the real 2x2 matrix X one can again obtain the SVD which leads to X = SxSxXqs- Since SxXcs — XqsSx we 
obtain the canonical decomposition $ = U 2 ° < i )CS ° f/i in terms of the symplectic transformations as stated in Eq. ( |2"2| ), with 

Xc-Xcs, Yc = Ybs, M 1 = S X , M 2 = ® Y S Y ® X , (49) 

demonstrating that (up to unitaries) the physical action of $ in this case corresponds to the action of <1> CS . Thus, we extended 
the theorem to lower rank cases of X and Y. □ 

Derivation of simplified expression for the one-shot Gaussian capacity 

In the following we show that the one-shot Gaussian capacity of a single-mode Gaussian channel $ can be expressed as 

C£($, N) = .max { X G (P, v) \ Tr[Mn + V moA ] <2N + 1}, (50) 



v in ,v m 



1\ / 1 



X G = 9[9-~j-9\v--y (51) 

g(x) — (x + 1) log 2 (a; + 1) - xlog 2 x, (52) 

where g(Q) = 0, Vj n is the CM of a pure Gaussian input state fulfilling det 2V m = 1, V mo d is the CM of a classical Gaussian 
distribution used to displace the input state and to generate the modulated input state with CM V[ n = V[ n + V mo< i where 
Tr[Vi n ] < 27V + 1. Furthermore, v = V det V ou t and D = \J dct V ou t are me symplectic eigenvalues of the output and 
modulated output state with CM V^, ut = $(VS n ) and V^ ut = $>(Vi n ), respectively. Equation ( f50| > states that among all possible 
Gaussian sources, where a source is characterized by a measure {j, G (da, dV p ) over the set of pure Gaussian states p G {a. V p ) 
of mean a an CM V p fulfilling det (2V P ) — 1, the source optimizing the Gaussian capacity of a given Gaussian channel 
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corresponds to using a single pure Gaussian state p G (0, V[n) with covariance V[ n fulfilling det (2Vj n ) = 1, modulated following 
a Gaussian bivariate distribution with CM Vmod- 

To achieve our goal we use the fact that the maximization inside the Gaussian capacity definition 



max 



5($[^ G ])~ / p G (da,dV p )S(*\p G (a,V p )]) 



(53) 



can be divided into two different steps. In the first step, among all the sources fi G (da, dV p ) belonging to the set T G G sharing 
the same average output state 



J (j, G (da,dV p )p G (a,V p ), 



we maximize the modified Holevo quantity 



X($,NJ G ) = S(<1>[P G })- min tl G (da,dV p )S(*\p G (a,V p )]) 



(54) 



(55) 



In the second and final step we optimize x($, N, p G ) over the average input state p G satisfying the energy constraint N, thus 
obtaining C G (&,N), It is well known that the minimum of the average output entropy appearing in equation (55 1 can be 

rewritten as the Gaussian entanglement of formation £' g [ctbe] of a given bipartite mixed state obe = U^pg (& |0) 



resulting from sending the average input mixed state p G through the unitary dilation [/$ of channel $, i.e. 

min f p G (da,dV p )S(^[p G (a,V p )}) = E G [a BE ], 



\ E ul 



(56) 



where the equality Tt^ctbe] = &[f> G ] is guaranteed by the Stinespring dilation of the channel. Indeed, the Gaussian entangle- 
ment of formation is defined as 



E G [a BE ] =E G [U^p JG ® |0) 



L/ll = min 



u G (da, dV p )E[U^p G (a, V p ) <g> |0) (0| U, 



(57) 



Here E[p G (a, V p )\ corresponds to the entanglement of a pure bipartite Gaussian state, which is uniquely quantified by the von 
Neumann entropy of any of its two reduced density operators. This connection, not only simplifies the capacity definition to 



max 



S(t>[p- G })-E G [UzP G ®\0) <0|t4] 



but also leads to the proof of Eq. (|50|>. In BOl it was shown that 



E G [p G (a,V)} = mm{E(p G (0, V P ))\V P < V}, 



(58) 



(59) 



where the infimum is taken over pure Gaussian states of CM V p . This implies the existence of a covariance matrix V p such that 
the output entropy of pure Gaussian states p G (ot, V p ) achieves the minimum in equation (56 1 for any a. Therefore, we chose 
the state p G (a, V p ) as the symbol states, i.e. VJ n = V p . Then modulating p G (a, VJ n ) according to a Gaussian distribution with 
covariance matrix V mod — V[ D — VJ n generates a source with average input state p G saturating the bound of equation (55 1. Thus, 
the one-shot Gaussian capacity for a general n-mode Gaussian channel <f> simplifies to 



C G ($, N) = max {S($(Vi„ + V mod )) - S(*(V m )) \ Tr[V in + V mod ] < 2nN + 1}. 



(60) 



v^„,v mod 

The entropy of an n-mode Gaussian state p G (a, V) can be calculated in terms of the n symplectic eigenvalues v t of V, namely 



S(p G ( a ,V)) = Y / 
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(61) 



where g{x) is defined in Eq. (f50|>. Thus, the one-shot Gaussian capacity of a general n-mode Gaussian channel simplifies to 



(- 1 

^ 1 v 



9 v. 



Tr[Mn + V^d] < 2nN +1 , 



(62) 



where Vi and v\ are the symplectic eigenvalues of the CM Vout 
Eq. (|62|) simplifies to Eq. (|50b. 



$(Vi n ) and CM V out — &(V ln ). For the one mode case 
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FIG. 4. Using the Bloch Messiah reduction any multimode pure Gaussian state can be generated from the vacuum and a set of single mode 
squeezers Si preceded and followed by linear multi-port interferometers ©i, ©2. (a) Reduction of the scheme for a collection of channels 
(<j> TH )® n and (b) Reduction for the fiducial channel (<I> F )® n (valid above an input energy threshold). 



Proof of Corollary|2] 

Corollary 2. For a single-mode Gaussian channel $ with parameters (r ^ 0,y > 0), there exists a fiducial channel $ F as 
defined in Theorem^such that 



C G ($, N) = C G (^ s) ,N) = g (\r\N + y cosh(2 s ) + _ g L + M_J^ 



if N>Ntto = \ fe 2|s| + j^sinh(2| S | 



(63) 



Proof. Corollary[l]can be straightforwardly extended to conclude that C G (<i>, N) — C G ($ F , N). This means that we only have 
to evaluate the Gaussian capacity of the fiducial channel $ F in order to find the Gaussian capacity for $. In the following we 
find C (<E> F , N) explicitly above the input energy threshold Nthr- 

The proof is structured as follows. First, we prove that the Gaussian minimum output entropy of thermal channels <I> TH is 
additive (corresponding to $ F with s = 0). Then we extend this proof to the fiducial channel for input energies above the input 
energy threshold 7V t hr (we present a simple and physically motivated proof, which is an alternative to the one in [ 17]). Then we 
show that this also implies the additivity of the one-shot Gaussian capacity in this energy domain. Finally, we derive the exact 
expression for the one-shot Gaussian capacity. 

In PTI it was shown that any pure n-mode Gaussian (input) state can be generated from a set of n vacuum modes, using 
n single-mode squeezers Si preceded and followed by a linear multi-port interferometer, corresponding to passive symplectic 
transformations ©1 and ©2 [see Fig.|4](a)]. The multi-mode vacuum state with CM 1/2 (where I is the 2n x 2n identity matrix) 
remains unchanged under the action of the first interferometer ©1 therefore, we can drop it. The action of the channel ^(y^) in 
terms of symplectic transformation then reads 

Kut = ^X m ® 2 SIS®]X Tli + r TH , (64) 

where S — ®™ =1 Si corresponds to the joint operation of single mode squeezers, ©2 is the symplectic transformation of the 
second linear multi-port interferometer and 



*TH = V M J \ FTH = yL (65) 



i=1 v sgnfrVK, . =1 
By inserting Xth^th = I between ©2 and S in Eq. d64} we obtain 



Kut = ©' QxthMSXth + y th) ®'\ (66) 
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where 0' = X T h©2^th and 0'Y TH Q' T = Y TK . One can confirm easily that 0T20' T = O and 0'0' T = I, so 0' is 
just another interferometer at the output. Thus, the general Gaussian input state entering the channel ($ TH )®" i s reduced to a 
product state [left hand side of Fig.|4](a)]. 

Now we calculate the entropy of the output state with CM V out - Since 0' does not effect its entropy we can drop it. Therefore, 
we are left with n vacuum modes entering the set of single-mode squeezers and then individually passing each channel $ TH . 
From this follows that the output entropy is additive and thus, 

min = min \ S ( (*™)°" fesis]) = min \ £ S(*™(V ini )) = min L - ±) , (67) 



v ln n v v ln n V \2 J J v ini n ^ v v v int n^ J \ 1 2 

where Vi n j = hSf is the CM of vacuum mode i exiting the single mode squeezer Si and Vi = y dct (<I> TH (Vi n j)) are the 
symplectic eigenvalues of the CM of the individual output states. The output entropy is minimized for V\ ai = 1/2, Vi, i.e. by 
fixing all one-mode squeezers S, = I. This leads to 

min -S(Kut) = 9 (y + ^P 1 ") ■ (68) 



----- ~ v ■ uui>/ C7 i a 1 

v^„ n \ 2 

Now let us consider the same problem for the fiducial channel <i> F as depicted in Fig. |4](b). We use again the Bloch-Messiah 
decomposition to decompose the general multi-mode Gaussian input state. The first interferometer ©i can again be dropped 
because it does not effect the n mode vacuum state. From the definition of the fiducial channel we have the equivalence 
Xp = Xth> = SX^hS, with S = diag(e s , e~ s ). This leads to the equality 

$ F (V) = SiXmS- 1 VS~ x Xra + Ytb)S. (69) 

As a consequence we can replace each fiducial channel by a thermal channel preceded by an anti-squeezer and followed by a 
squeezer [see right hand side of Fig.|4](b)]. 

Now we focus again on the minimization of the output entropy. Then, the set of squeezers S can be dropped since they do 
not change the entropy and we have no energy constraint on the output. We showed above that the entropy for the joint map 
(<& TH )®™ is minimized by the n-mode vacuum state. Thus, the multi-mode input state that minimizes the output entropy of 
the fiducial channel has to be in the n-mode vacuum state after passing the n anti-squeezers S^ 1 [see right hand side of Fig. [4] 
(b)] . Therefore, one fixes the interferometer ©2 = I and choses each squeezer Si at the input to undo the anti-squeezers, i.e. 
Si = S, Thus, the entropy minimizing n-mode Gaussian input state for the channel ($ TH )® n is the product state with CM 
V in = ®f =1 S 2 /2. It follows that 

min '-S ((^f n (V in )) = 9 (y+ M^i) , (70) 

where the minimization requires a certain amount of energy to undo each squeezer, which will be taken into account later. 

We show now that Eq. ( f70] i leads to the additivity of the one-shot Gaussian capacity of <f> F above the input energy threshold 
iVthr. Using the expression for the one-shot Gaussian capacity C^($ F , N) stated in Eq. (|60j> and using the definition of the 
Gaussian capacity 

C G ($,A0= lim ~C9(9 9n ,nN), (71) 

we can state the following upper bound: 

C G ($ F , N) < max S ($ F (M„ + V mod )) - lim min -S ( ($ F f "(Mn)) , (72) 

where the first term only needs to be maximized for a single use of the channel due to the subadditivity of the entropy. It is 
known that a thermal state maximizes the entropy, therefore, the optimal modulated output state is a thermal state carrying the 
total number of photons, i.e. <I> F (Vi n + V, noc i) = V^ ut = diag(i^, P), where 

v = |r|iV + ?/cosh(2s) + \t\/2. (73) 

The second term in Eq. f72] > was already evaluated above [see Eq. ((70])]. In summary, we found 

max 5($ F (Mn + K nod )) = g (\t\N + ycosh(2 S ) + ^-^-) , (74) 

Vin.Vmod \ Z y 

lim rmn±s((**f n (V in ))=g(y + ^). (75) 
rn-oo vin n V /V 2 / 
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The encoding which realizes both the maximum and the minimum in Eqs. |74} a ( |75j ) is given by 

{N + h-& sinh(2s) \ 

V in = S 2 /2, , Knod = Mn - Mn, M„ = 2 ^ ^ x „ .,,„,■ (76) 

\ N + ^ + \^\ smh(2s) J 



This encoding can only be realized if 



N > N thr = \ (e 2 ^ + sinh(2| S |) - 1 j . (77 ) 



because for N < Nthx it would imply Vmod < which would be non-physical. Thus, we have shown that 

C G ($ F ,iV) = C«($ F ,A0, N>N thl , (78) 
where C^($ F , N) is the right hand side of Eq. ( |63"j ). This proves the Corollary. □ 

Additional upper bounds on the classical capacity 

The upper bound on the classical capacity stated in Corollary [3] was obtained by generalizing the bounds that were found for 
thermal channels $ TH (with r > 0) in ll24ll . Recently, additional upper bounds were obtained for the same channels |39| and we 
extend them now to general channels $ with r > 0, y > 0. 

The bounds were obtained by maximizing the first term of the Holevo quantity [defined in Eq. ( fT0| ) in the main text] and by 
obtaining a lower bound b on the second term, i.e. 

lim -min / fi(dx) S (($ TH ) S) ™[/5 2; ]) > b, (79) 

n— >oo n fJ. J 

where in total six bounds b are stated in (39). We stated in the proof of Corollary [2] that the fiducial channel is equivalent to 
a thermal channel preceded by an anti-squeezer and followed by a squeezer [see Fig. |4](b)]. The following squeezer does not 
change the output entropy. Furthermore, one can always undo the preceding squeezer because the bound b is not subject to an 
energy constraint. Therefore, any bound b on the minimal output entropy of the thermal channel is as well an upper bound on the 
minimal output entropy of the fiducial channel. In addition, the first term of the classical capacity is known to be maximized by a 
one-mode Gaussian thermal state carrying the total number of photons of the system. This quantity is already given in Eq. ( |74] i, 
i.e. 

lim - max S (($ F )®"[^]) = g (\t\N + ycosh(2s) + J • (80) 

Corollary [j] states that C($, N) = C($ F , N). Therefore, any bound on the classical capacity of the fiducial channel $ F is also 
a bound on the classical capacity of a an arbitrary channel $. Thus, we obtain a list of upper bounds on the classical capacity 
that reads 

C%N) <g^\T\N + yco S h(2 S ) + ^^j - b, r > 0,y > 0, (81) 

where b has to be taken from [39] and substituted in the latter. Note that with increasing s those bounds become less and less 
tight because b does not depend on s. However, as in the case of thermal channels lf39l some of those bounds in a certain region 
of channel parameters are tighter than the bound we presented in Corollary [3] 



